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Background: Multiple-locus variable-number tandem-repeat fingerprinting (MLVF) is based 
on multiplex PCR, utilizing variable number tandem repeat. Our goal was to compare the 
performance of MLVF in distinguishing clinical Staphylococcus aureus isolates with that of 
pulsed-field gel electrophoresis (PFGE), which has traditionally been the gold standard. 

Methods: Sixty-three clinically significant S. aureus isolates were tested using both PFGE 
and MLVF. Multiplex PCR for MLVF was performed using PCR primers for clfA, clfB, 
sdrCDE, sspA, and spa. PFGE was performed with genomic DNA fragments generated by 
Sma\ endonuclease digestion. Banding patterns of MLVF or PFGE were analyzed using 
InfoQuestFP software. 

Results: The hands-on time of our modified method was about 3 h, on average, for each 
of 18 isolates. PFGE (80% cutoff) or MLVF (75% cutoff) separated all of the 63 isolates into 
13 and 12 types, respectively. Three types generated by PFGE were identical to those gen- 
erated by MLVF. PFGE and MLVF yielded similar Simpson's diversity indices, indicating 
similar discriminatory power. The overall concordance between PFGE and MLVF was low, 
as represented by adjusted Rand indices (0.266-0.278). PFGE predicted MLVF type better 
than MLVF predicted PFGE type, as reflected by Wallace coefficients (PFGE cutoff 80% vs. 
MLVF cutoff 75%, 0.389 vs. 0.233). Analysis of the relationship between a pair of isolates 
showed 91.0% concordance between the PFGE (80% cutoff) and MLVF (75% cutoff). 

Conclusions: Our simple, low-cost, modified MLVF protocol can effectively discriminate 
between S. aureus clinical isolates. MLVF can replace PFGE for the hospital infection con- 
trol of S. aureus. 

Key Words: Staphylococcus aureus, Pulsed-field gel electrophoresis, Multiple-locus vari- 
able-number tandem-repeat fingerprinting 
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INTRODUCTION colonized with Staphylococcus aureus, colonization increases 

the risk of S. aureus infection when the host immune system is 
Although 20-30% of the general population is asymptomatically compromised [1]. S. aureus infections vary from superficial 
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wound infections to invasive diseases, such as deep abscesses, 
osteomyelitis, and bacteremia [2]. Methicillin-resistantS. aureus 
(MRSA) is of especially great concern because MRSA infection 
extends the length of hospital stay and increases antibiotic use, 
costs, and mortality [3]. 

Traditionally, MRSA strain typing was accomplished by pulsed- 
field gel electrophoresis (PFGE) [4], However, the operational 
time of PFGE is at least 72 hr, the cost is relatively high, and the 
technique requires specialized training [5]. Sabat et al. [6] there- 
fore designed a multiplex PCR scheme called multiple-locus 
variable-number tandem-repeat fingerprinting (MLVF), using 5 
variable number tandem repeat (VNTR) loci. Studies have shown 
that MLVF can distinguish among MRSA strains and simultane- 
ously provide results that parallel those of PFGE [5, 7, 8]. 

We aimed to investigate the usefulness of MLVF for typing S. 
aureus isolates of clinical significance. We also tried to modify 
the MLVF method described previously, in order to maximize 
the technical advantage of MLVF, in terms of time, cost, and 
simplicity, compared with PFGE. 

METHODS 

1. Bacterial isolates 

Sixty-three hospital-acquired S. aureus isolates recovered in 
2009 were selected according to Centers for Disease Control 
and Prevention (CDQ/National Healthcare Safety Network 
(NHSN) criteria [9]. Isolates that were regarded as contaminant 
or non-pathogenic were excluded. Among the 63 isolates, 15 
were methicillin-susceptible, and 48 were methicillin-resistant. 
Most of the bacterial strains were isolated from blood culture 
and groups of 2 or 3 isolates that were considered epidemiologi- 
cally related were chosen for analysis. The isolates were recov- 
ered from the following clinical sources: blood culture (N=43), 
ascitic fluid (N = 2), cerebrospinal fluid (N=3), sputum (N=3), 
endotracheal aspiration (N=3), pus (N=2), T-cannula tip (N=2), 
wounds (N =2), throat (N = 1), tissue (N = 1), and drain (N = 1). 

2. DNA preparation 

Bacterial isolates were subcultured in 5% sheep blood agar over- 
night at 37°C. Total genomic DNA was extracted from S. aureus 
isolates using the QIAamp DNA Mini Kit (Qiagen, Hilden, Ger- 
many), according to the manufacturer's instructions. Purified 
DNA was stored at -20°C. 

3. MLVF typing 

The multiplex PCR assay containing clfA, clfB, sdrCDE, spa, and 



sspA primers was performed as previously described by Sabat 
et al. [6], with modifications. For convenience, the AccuPower 
HotStart PCR PreMix kit (Bioneer Co., Daejeon, Korea) was 
used in the PCR reaction. The kit contained a mixture of 1 unit 
of HotStart DNA polymerase, lx PCR buffer, 250 uM of each 
dNTP, and 1.5 mM MgCI 2 in each reaction tube. A mixture of 
the following concentrations of each primer and 1 uL of tem- 
plate DNA (20-50 ng) was added to the reaction tube: 0.15 uM 
of clfA-F (forward) and clfA-R (reverse), 0.2 uM of clfB-F and 
clfB-R, 0.15 uM of sdrCDE-F and sdrCDE-R, 0.1 uM of spa-F 
and spa-R, and 0.3 uM of sspA-F and sspA-R. The thermal cy- 
cling was performed in a PTC-100 Thermal Cycler (MJ Research, 
Waltham, MA, USA) as follows: 94°C for 5 min, followed by 30 
cycles at 94°C for 30 sec, 55°C for 30 sec, 72°C for 1 min, with a 
final extension at 72°C for 7 min. 

4. Visualization of PCR products 

The amplified PCR products and 100-bp DNA ladder marker 
(Roche Applied Science, Indianapolis, IN, USA) were resolved 
by electrophoresis in a 2% agarose gel in 0.5 x Tris-borate-EDTA 
(TBE) buffer at 50 V for 80 min using Mupid-2plus (Takara Bio 
Inc., Otsu, Japan). The gel image was saved as a TIFF file and 
analyzed using InfoQuestFP Software (Bio-Rad Laboratories 
Inc., Hercules, CA, USA) under the following tolerance settings: 
optimization, 0.5%, and position tolerance, 1.25% [5]. The band 
patterns were clustered into MLVF types using Dice's coefficient 
and the unweighted pair-group method with arithmetic mean 
(UPGMA). The effect of grouping similar MLVF patterns together 
into a smaller number of MLVF types was examined using cut- 
offs of >75% and >70% relatedness. 

5. PFGE 

Genomic DNA of the isolates was prepared, digested by restric- 
tion enzyme Sma\ (Roche Diagnostics GmbH, Mannheim, Ger- 
many), and separated in a GenePath system (Bio-Rad Labora- 
tories Inc.) according to the CDC PulseNet protocol with modifi- 
cation [10]. Sma\ PFGE patterns were saved in TIFF files. The 
file was exported to the InfoQuestFP software (Bio-Rad Labora- 
tories Inc.) and analyzed using the Dice coefficient-UPGMA. A 
dendrogram was generated to examine the relatedness of PFGE 
patterns for all study isolates, and cutoff levels of 80% and 75% 
were applied to this dendrogram. 

6. Calculation of concordance 

Simpson's index of diversity was calculated to measure the dis- 
criminatory power of the typing systems. This index indicates 
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the probability that 2 strains sampled randomly from a popula- 
tion will belong to 2 different types [11]. The formulas for Simp- 
son's index (D) and the confidence interval (CI) are presented 
in the following equations: 



</-{-[i*t-<Li!] 

CI = [D-2V?, D+2V0 2 ] 

where N is the total number of strains in the sample population, 
S is the total number of types described, m is the number of 
strains belonging to the jth type, and n, is the frequency n/N. 

To compare the 2 sets of results produced by PFGE and MLVF, 
we used 2 clustering comparison coefficients (adjusted Rand 
and Wallace) generated by a script from http://biomath.itqb.unl. 
pt/ClusterComp, as previously described [11]. We also evaluated 
the concordance of each pair of isolates as previously described 
[12]. All possible pairs of isolates were cross-classified on the 
basis of matched or mismatched types. The resulting 2x2 table 
was evaluated using the chi-square statistic, and the percentage 
of concordant cells was calculated. 

RESULTS 

1. Technical aspects of MLVF 

Using a commercial kit containing a premixture of PCR reagents, 
we were able to perform MLVF more easily by simply transfer- 
ring the template DNA and primers into the reaction tubes. This 
MLVF protocol, modified from a previous study [6], is sufficiently 
simple for use in clinical laboratories, and 18 isolates can be 
typed in a single working day if pure cultures of bacterial iso- 
lates are available. The hands-on time for 18 tests was approxi- 
mately 3 hr on average (2 hr for DNA extraction, 30 min for 
preparation of PCR reaction, and another 30 min for MLVF pat- 
tern reading), with most of the time spent on the extraction of 
DNA. Experienced personnel would be able to analyze the 
MLVF patterns by computer software within 30 min. 

2. Molecular typing by PFGE and MLVF 

All 63 S. aureus clinical isolates could be typed by PFGE (Fig. 1) 
and MLVF (Fig. 2). Because defining clonality in the PFGE using 
a similarity value of 80% as a cutoff is considered the gold stan- 
dard [13], we used 80% and 75% cutoff values for PFGE and 



75% and 70% for MLVF. We also tried an 80% cutoff for MLVF; 
however, we found that some isolates were not grouped into the 
same type predicted by PFGE (data not shown). PFGE distin- 
guished the isolates into 13 or 14 types and an additional 24 or 
20 unique patterns at cutoff values of 80% and 75%, respectively 
(Fig. 1). Similarly, MLVF separated the 63 S. aureus isolates into 
12 or 13 types and additional 16 or 12 unique patterns at cutoff 
values of 75% and 70%, respectively (Fig. 2). PFGE types P10 
and P13 were identical to MLVF types Ml and M5, respectively. 
A very good correlation was observed between PFGE types P8, 
P5, and P6 and MLVF types M3, M4, and M9, respectively. In 
these types, only 1 or 2 isolates composing each MLVF type 
were different from the isolates composing the PFGE type. How- 
ever, MLVF type M8, which was composed of 9 isolates at both 
75% and 70% cutoff values, was divided into 5 different PFGE 
types, indicating that MLVF had a poorer discriminatory power 
than did PFGE. 

3. Comparison of typing methods 

To compare the discriminatory power of PFGE and MLVF, we 
determined the mean number of isolates per type and Simpson' 
s diversity index for all the isolates (Table 1). MLVF (75% cutoff) 
and PFGE (80% cutoff) yielded similar D-values of 0.954 (95% 
CI 0.933-0.974) and 0.972 (95% CI 0.957-0.988), respectively. 

To compare the congruence between type assignments using 
PFGE and MLVF, we calculated adjusted Rand indices (Table 2) 
and Wallace coefficients (Table 3) for each cutoff value of PFGE 
and MLVF. The overall congruence between PFGE and MLVF 
was presented by adjusted Rand indices, ranging from 0.266 to 
0.278. As expected from the Wallace coefficients, PFGE pre- 
dicted MLVF type better than MLVF predicted PFGE type (Wal- 
lace coefficients, PFGE cutoff 80% vs. MLVF cutoff 75%, 0.389 
vs. 0.233; PFGE cutoff 80% vs. MLVF cutoff 70%, 0.426 vs. 
0.232; PFGE cutoff 75% vs. MLVF cutoff 75%, 0.359 vs. 0.256; 
PFGE cutoff 75% vs. MLVF cutoff 70%, 0.391 vs. 0.253). 

Although the congruence between PFGE and MLVF in typing 
these clinical isolates was poor, notable agreement was found 
when the relationship between a pair of isolates was analyzed 
(Table 4). Cross-classification of the isolates, based on matched 
or mismatched schemes by PFGE (80% cutoff) and MLVF (75% 
cutoff), showed that the 2 typing systems were 91.0% concor- 
dant. Using different cutoffs, the PFGE and MLVF showed 88.3% 
to 90.1% concordance (data not shown). Any given pair of iso- 
lates distinguished by one method tended to be distinguished 
by the other. 
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Fig. 1. Pulsed-field gel electrophoresis (PFGE) types of the study isolates generated by the unweighted pair-group method with arithmetic 
mean (UPGMA) algorithm, using 80% and 75% similarity cutoff values. The isolates that were not included in the above types were con- 
sidered to be unique patterns. Isolate numbers presented in boldface were not included in any types generated by the 80% cutoff. 



DISCUSSION 

Several studies based on VNTR schemes have been published 
forS. aureus strain typing [5-8, 14-16]. Most of the studies, how- 
ever, were performed using selected MRSA isolates of certain 
limited PFGE types, which might introduce a selection bias. In 
this study, we chose clinically significant S. aureus isolates to 



evaluate the utility of MLVF, with respect to infection control in a 
hospital setting. Because PFGE is labor intensive, time consum- 
ing, and expensive, routine use of this procedure in a clinical 
laboratory is rather demanding. Therefore, we tried to determine 
whether PFGE could be replaced with the more "simplified" 
MLVF. We used a commercial premix PCR kit, instead of pre- 
paring a mixture of PCR reagents before amplification. After PCR 
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Fig. 2. Multiple-locus variable-number tandem-repeat fingerprinting (MLVF) types of the study isolates generated by the unweighted pair- 
group method with arithmetic mean (UPGMA) algorithm using 75% and 70% similarity cutoff values. The isolates that were not included in 
the above types were considered to be unique patterns. Isolate numbers presented in boldface were not included in any types generated 
by the 75% cutoff. 



amplification, we performed electrophoresis for only 80 min with 
a 12-cm length system, which can easily be utilized in clinical 
laboratories. Our modified method could save much time, com- 
pared to PFGE, which requires approximately 3 days (2 days for 
DNA agarose disc preparation, followed by 1 day of PFGE), as 
well as the previous MLVF procedure that utilized a master mix, 



which was prepared by adding each ingredient to each reaction 
tube, plus electrophoresis using an extended gel, which re- 
quired 3 hr 30 min for running [6]. Even the cost of our modified 
MLVF was much lower than that of PFGE ($10 vs. $15 per test). 

As shown by Simpson's diversity index, in our study, MLVF 
showed discriminatory power similar to that of PFGE. This find- 
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Table 1. Simpson's diversity indices of the genotyping methods for Table 2. Adjusted Rand index for PFGE and MLVF using different 



all isolates typed 



Methods N 


o. of types 


No. of isolates 
per type (mean±SD) 


D (95% CI) 


PFGE (80% cutoff) 


13 


3.0 ± 1.47 


0.972 (0.957-0.988) 


PFGE (75% cutoff) 


14 


3.1 ± 1.44 


0.967 (0.951-0.983) 


MLVF (75% cutoff) 


12 


3.9 ± 1.98 


0.954 (0.933-0.974) 


MLVF (70% cutoff) 


13 


3.9 ± 2.02 


0.949 (0.929-0.970) 



Abbreviations: D, Simpson's index; CI, confidence interval; PFGE, pulsed- 
field gel electrophoresis; MLVF, multiple-locus variable-number tandem-re- 
peat fingerprinting. 



cutoffs 



Typing method 

PFGE cutoff 80% 
PFGE cutoff 75% 
MLVF cutoff 75% 
MLVF cutoff 70% 



Adjusted Rand index 



PFGE cutoff 80% PFGE cutoff 75% MLVF cutoff 75% 

1 

0.860 1 

0.266 0.271 1 

0.275 0.278 0.950 



Abbreviations: PFGE, pulsed-field gel electrophoresis; MLVF, multiple-locus 
variable-number tandem-repeat fingerprinting. 



Table 3. Wallace coefficients for PFGE and MLVF using different cutoffs 



Typing method 




Wallace coefficient (95% CI) 




PFGE cutoff 80% 


PFGE cutoff 75% 


MLVF cutoff 75% 


MLVF cutoff 70% 


PFGE cutoff 80% 




0.944 (0.886-1.000) 


0.389 (0.275-0.503) 


0.426 (0.308-0.544) 


PFGE cutoff 75% 


0.797 (0.706-0.888) 




0.359 (0.254-0.465) 


0.391 (0.276-0.506) 


MLVF cutoff 75% 


0.233 (0.145-0.321) 


0.256 (0.161-0.350) 




1.000 (1.000-1.000) 


MLVF cutoff 70% 


0.232 (0.146-0.319) 


0.253 (0.158-0.347) 


0.909 (0.836-0.982) 




Abbreviations: CI, confidence interval; PFGE, pulsed-field 


gel electrophoresis; MLVF, multiple-locus variable-number tandem-repeat fingerprinting. 



Table 4. Cross-classification of all possible pairs of isolates based 
on matched or mismatched PFGE and MLVF types 



PFGE 

(80% cutoff) 


MLVF (75% cutoff) 


Sum 


Match 


Mismatch 


Match 


31 


30 


61 


Mismatch 


145 


1,747 


1,892 


Sum 


176 


1,777 


1,953 



P< O.OOl. 

Concordance = 3 |t, 1 J^ 47 xl00=91.0(%) 
1,953 

Abbreviations: PFGE, pulsed-field gel electrophoresis; MLVF, multiple-locus 
variable-number tandem-repeat fingerprinting. 



ing is consistent with the findings of a previous study that used 
8 different VNTR loci, including spa and sspa [14]. However, the 
congruence between types defined by PFGE and MLVF was 
fairly low, as reflected by adjusted Rand index and Wallace co- 
efficients. The adjusted Rand index allows a quantitative evalua- 
tion of the global congruence between the 2 partitions [11]. On 
the other hand, Wallace coefficients, given a particular typing 
method, provide an estimate of how much new information is 
obtained from another method. The adjusted Rand index be- 
tween PFGE and MLVF calculated in our study was similar to 
that calculated in 2 previous studies [16, 17]. However, the Wal- 
lace coefficient of PFGE predicting MLVF type in our study was 
lower than that calculated in the 2 previous studies. The dis- 



crepancy in the results may be explained by the difference in 
the S. aureus isolates used. Since the electrophoretic pattern of 
MLVF is based on staphylococcal interspersed repeats, the dif- 
ference in the bacterial isolates might have affected the cluster- 
ing profile. Those previous studies enrolled EMRSA-15 (ST22) 
and EMRSA-16 (ST36/ST30) or ST398 MRSA. In our study, we 
used various clinical isolates with diverse PFGE or MLVF pat- 
terns that could be classified into 37 types by PFGE (80% cut- 
off) and into 28 types, including unique types by MLVF (75% 
cutoff). Moreover, 15 of 65 isolates were methicillin-susceptible 
S. aureus (MSSA), which is known to be more polyclonal than 
MRSA [18]. 

Another explanation for the relatively low Wallace coefficients 
of PFGE predicting MLVF types in our study could be that we 
used only Sma\ for DNA digestion in the PFGE. In the study by 
Rasschaert et al., BsfZI, SacW, Apa\, and combinations of these 
enzymes were used [17]. The combination of all 3 restriction 
enzymes yielded a Wallace coefficient of 1.00 in predicting MLVF 
types in that study. 

However, when the relationship between 2 specific isolates 
was analyzed, PFGE (cutoff 80%) and MLVF (cutoff 75%) showed 
91.0% concordance, calculated by each cell (Table 4). In other 
words, every pair of isolates can be distinguished by either MLVF 
or PFGE. Although the InfoQuestFP software provided excellent 
dendrograms by applying multilinkage where branches degen- 
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erate, every dendrogram would confront loss of information while 
converting the data set of similarity values into a simplified hier- 
archical structure [19]. Furthermore, similar to PFGE analysis, 
MLVF also required careful visual inspection to confirm the re- 
sults owing to small band shifts. 

Even though the Wallace coefficients for MLVF predicting 
PFGE type (gold standard cutoff of 80%) were similar at both 
MLVF cutoffs of 75% and 70% (0.233 vs. 0.232), we suggest that 
a 75% cutoff be used for MLVF, since it showed higher discrimi- 
natory power. 

In conclusion, our modified protocol of MLVF demonstrated 
good performance in typing of clinically significant S. aureus 
isolates selected on the basis of the CDC/NHSN criteria. Further- 
more, we were able to shorten the running time by using a PCR 
premixture kit and performing electrophoresis using a 12-cm 
length system. Because MLVF has large advantage over PFGE 
in terms of technical aspects, MLVF could be used in place of 
PFGE, especially when urgent investigation of a hospital S. au- 
reus outbreak is required. 
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